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The fluctuation- dissipation relation is calculated for a class of stochastic models obeying 
a master equation. The transition rates are assumed to obey detailed balance also in the 
presence of a field. It is shown that in general the linear response cannot be expressed via 
time-derivatives of the correlation function alone, but an additional function £(t,t w ) , which 
has been rarely discussed before is required. This function depends on the two times also 
relevant for the response and the correlation and vanishes under equilibrium conditions. It 
can be expressed in terms of the propagators and the transition rates of the master equation 
but it is not related to any physical observable in an obvious way. Instead, it is determined 
by inhomogeneities in the temporal evolution of the distribution function of the stochastic 
variable under consideration. £(t,t w ) is considered for some examples of stochastic models, 
some of which exhibit true non- equilibrium dynamics and others approach equilibrium in 
the long term. In particular, models in which a relevant variable, e.g. a magnetization, is 
related in a prescribed way to the states of the system are considered as projections from 
a composite Markov process. From these model calculations, it is conjectured that £(t,t w ) 
vanishes when one is concerned with measurements of the analogue of a structure factor 
for large wave-vectors or if every transition among the states randomizes the value of the 
stochastic variable considered. 

PACS Numbers: 64.70 Pf,05.40.+j,61.20.Lc 



I. Introduction 

The out-of-equilibrium dynamics of stochastic models has gained intensive interest in the 
last decade. In particular, the deviations from the fluctuation dissipation theorem (FDT), 
relating the linear response to the two-time correlation function have been investigated 
in great detail, for a recent review seepQ. Starting with an investigation of the spherical 
p-spin-glas modelp], much attention has been paid to study the behavior of the response 
and the correlation for models of glassy dynamics. While in equilibrium the FDT relates 
the response to the correlation in a unique way, this does not hold in out-of-equilibrium 
situations. The violations of the FDT usually are parameterized via the introduction of a 
function X(t,t w ), which is defined via: 
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Here, the correlation function of a quantity M(t) is defined by C(t, t w ) = (M(t)M(t w )) and 
the corresponding response to a field conjugate to M{t) is R(t,t w ) = 5(M(t)) /5H(t w )\ H=0 
for t>t w . In case that M is a so called neutral variable"""), X(t,t w ) is independent of M 
and the long time limit 1^ allows the definition of an effective temperature ["£]. 

The value of has been calculated for a variety of models with different results. One 
class of models that have been considered are coarsening models [5 , for which is known 
to vanish"""). Examples of such models are the well known spherical modelpH IB] and the 
O(N) model in the limit of large iV[9 3 , as well as the Ising models in one [Til] or higher 
dimension In the context of models for glassy dynamics for some discontinuous mean 
field spin models a different behavior has been found"!""]. in particular, a relation between 
the degree of replica symmetry breaking and X^ has been established. In addition, in some 
time sectors, in which both, the response and the correlation obey some scaling relations, 
it has been found that X(t, t w ) is a function of the correlation alone, X(C). For disordered 
systems, a connection between the out-of-equilibrium dynamics and static properties has 
been established [I""]. For coarsening systems, however, such a relation does not exist """""] . 
In addition to these analytical calculations, a number of molecular dynamics simulations 
have been performed on model glassforming liquids, for a recent review see ref.|15j. 

Some of the quoted models are soft-spin models in the sense that the stochastic dy- 
namics is calculated from a Langevin equation. The classical treatment of FDT violations 
for unfrustrated stochastic models with Langevin dynamics has been given in ref.|16j. 
where various examples, including a simple random walk have been considered. In addi- 
tion to models obeying a Langevin equation, the out-of-equilibrium dynamics of models 
with a dynamics determined by a master equation [T7j have been investigated, in particular 
in the context of the aging dynamics in spin glasses. In this context also the function 
£(t,t w ), which will be the central topic of the present paper, has been discussed for the 
first time |18l I19j . A well known model in this context is Bouchauds trap model|20] and a 
number of investigations of the FDT violations have been presented |21 [ I22 [ \'2'A\. 

In the present paper, I consider models for which the probability distributions obey a 
master equation (ME) and I assume that the transition probabilities fulfill the conditions 
of detailed balance also in the presence of a perturbing field. The behavior of the function 

t w ) will be discussed in detail for some specific models. I will mainly focus on models, in 
which the dynamics of some stochastic variable M(t) is determined via transitions among 
the states of the system under consideration. Therefore, one has to assign values Mj. to 
the states k. This can be done in a variety of different ways. Here, I will discuss a class 
of models in which the stochastic process M(t) is viewed as a projection from a two- 
dimensional composite Markov process (CMP) [""""]. When the transition rates are chosen 
in an appropriate way, this procedure allows to treat several different models on the same 
footing, including random magnetic models. The outline of the paper is the following. 
In the next section the general formalism will be discussed. This includes a discussion of 
trap models because one generally finds that £(t,t w ) vanishes for these models [2"""] "2"2"] 12*3] . 
Section III is devoted to random magnetic models and a particular example of a kinetic 
random energy model is discussed in detail for illustrative purposes. Models for tagged 
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particle motion as they have been applied to the dynamics in supercooled liquids are 
considered in Section IV and the conclusions are presented in Section V. 

II. Master equations and FDT violations 

A. General formalism 

Throughout this paper a stochastic dynamics according to a master equation (ME)[U| is 
assumed. In a discrete notation let G k i(t) be the conditional probability to find the system 
in 'state' k at time t provided it was in 'state' I at time t = 0. At this point it is not 
necessary to specify the meaning of the term 'states'. Denoting the transition probabilites 
for a transition from state k to state I by Wi k , the ME reads as: 

Gki(t) = - E W nk G M (t) + W kn G nl {t) (2) 

n n 

Of course, the same ME is obeyed by the populations of state k, Pk(t) = J2i G k i(t)pi(0). In 
addition, I only consider transition probabilities that obey detailed balance 

WuV? = W lk p e k q (3) 

where the p e k q are the populations in thermal equilibrium. These are invariant with respect 
to the master operator, i.e. p e k q = Y^iGki(t)pi q - 

Throughout the present paper it is assumed that the system is prepared in some initial 
state described by a fixed set of populations, p k . The populations then evolve according 
to Pkif) = J2i Gkii^P^- In the general case, of course, the p k are different from the p e k , but 
they still fulfill the sum rule J2k Pk = ^ as required for probabilities. 

Central to the topic of the present paper is the two-time correlation function 

C{t, t w ) = (M(t)M{t w )} = MkMiGkiit - t w ) Pl (t w ) (4) 

k,i 

where M k is the value of M(t) in state 'fc'. In this expression, t w denotes the time that has 
evolved after the initial preparation of the system in the populations p k . In the following t > 
t w will always be assumed. In addition, the correlation function can always be decomposed 
according to 

C(t,t w ) = C eq (t-t w ) + AC(t,t w ) with C eq (t-t w ) = Y,M k M l G kl (t-t w )p e l q (5) 

k,l 

because of the properties of G k i(t w ). Note that this property only holds if the W k i obey 
detailed balance, cf. eq.©. From G ' k i(t w — >• oo) =p k q one sees that AC(t,t w ) vanishes for 
long waiting times, because exactly this long time limit has been subtracted from AC(t, t w ). 
Eq.© holds for all models considered in the present paper. 

If (M(t))^0, it is advantageous to consider C(t, t w ) = (M(t)M(t w )) - (M(t))(M(t w )), 
given by C(t, t w ) =£ M M k M x [G kl {t-t w ) - p k (t)} Pl (t w ). 
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In order to calculate the linear response of the system 

S(M(t)) 



R[t, 



5H(t v 



(6) 

H=0 



the dependence of the transition probabilities on a field H conjugate to M has to be 
fixed. In principle there is no restriction regarding this dependence. From equilibrium 
considerations one expects a Boltzmann-like dependence, e l3HM . This, however, does not 
fix the dependence on the values of M k in the initial or final state of a k — > I transition. In 
the present paper, I choose the following form, which assures that the system also in the 
presence of the field fulfills detailed balance: 

W kl (H) = W kl e /3HX ^ with X kl = aM k - (1 - a)M l (7) 

Here, a is a parameter that can take on any value. The same dependence has also been 
used by Bouchaud and Dean[2T] in a study of the aging properties of the trap model 20 . 
Very recently, Ritort has generalized this dependence in using X k i = ^M k — \iMi with 
arbitrary 7 and /i with the effect that the W k i(H) do no longer obey detailed balance[22 . 

Even though a in principle can take on any value, there often will be some guiding 
principle. For example, if the states denote the energies in a canonical ensemble, one 
expects that a can be determined from the dependence of the unperturbed W k i on k and 
If the W k i are of a form allowing a Kramers-Moyal expansion |24j. and therefore the ME 
has a well defined Fokker-Planck equation as a limit one would naturally choose a = 1/2. 
However, it has to be pointed out that usually the states k are understood as metastable 
states or components |25j in connection with glassy systems. Then the corresponding free 
energies are to be viewed as coarse-grained quantities [1^ and one does no longer have a 
strict relation of a to the unperturbed W k \. 

Using eq.(J7J), the linear response R(t,t w ), eq.©, can be calculated with the result 

R(t, t w )=pY, M k [G kn {t - t w ) - G kl (t - t w )} W nl X nlPl (t w ) (8) 

k,l,n 



■ dC(t,t w ) , dC(t,t w ) 

a — — (1 - a) — + a£{t,t v 

ot,„ at 



(9) 



After some algebra this can be related to the correlation function: 

R(t,t w )=(3 

where I defined the function 

f (*, *») = E M k MiG kn {t - t w ) [W nlPl (t w ) - W lnPn (t w )} (10) 

k,l,n 

It is important to point out that eq.(jnj) holds for arbitrary Markov processes obeying 
detailed balance. Additionally, the function £(t,t w ) cannot be related to a time derivative 
of the correlation function and therefore the response is not determined by C(t,t w ) alone 
in the general case. £(t,t w ) plays a similar role as the asymmetry in the treatment of 
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the response derived from a Langevin equation [TB]. Of course, a relation between R(t,t w ) 
and C(t,t w ) mentioned above is easily obtained from eq.Q by noting that dC(t,t w )/dt = 
dC(t,t w )/dt - d((M(t))(M(t w )))/dt and similarly for dC(t,t w )/dt w . 

If the system is prepared in an equilibrium state initially, pl=p e k 9 , one has Pi(t w ) = p^ q 
and eq.© shows that £ e <j(Mto) = 0- Furthermore, the response and the correlation are 
functions of the time-difference only, i.e. C eg (t,t w ) =C eq (t — t w ) and R eq (t,t w )=R eg (t — t w ) 
and thus 

ReM = -P^P- (11) 

which is just the well known FDT. 

To the best of the authors knowledge, eq.© has not been derived in this form before. 
However, equations similar to eq.(jHJ) have been given for various models in the literature. 
Hoffmann and Sibani^Hj have derived eq.Q for the special case a = 1. Later on Bouchaud 
and Dean[2I] give a similar expression with, however, £,(t,t w ) = 0. Furthermore, Fielding 
and Sollich[3] derived eq.® for the special case of a = 0. In all these cases the authors 
considered model systems where the quantities denote a magnetization which is at- 
tached to state k in some way. Such models will be treated as special cases in the following 
chapters. Before doing so it is instructive to have a somewhat closer look at the quantity 
£(t,t w ). As noted by Hoffmann and Sibani, £(t,t w ) vanishes if the relaxation to equilib- 
rium is determined by distribution functions that are equilibrated with respect to the states 
k and depend on time only parametrically. In order to see this explicitly let us assume 
that the initial conditions are such that one can write Pk{t w ) = p1 9 5k{t w ) with unspecified 
functions 6k(t w ). In this case eq.fITU]) reads as 

i(t,t w ) = ]T MkMtGUt ~ t w )W nlV T [SiiU) - 5 n (t w )] 

From this expression it is evident immediately that £(t,t w ) vanishes for Sk{t w )=5(t w ) V7c. 

In the following sections some specific choices of the states k and of the dynamic vari- 
ables will be presented and eq.© will be discussed for these cases. In particular, I will 
concentrate on the function £(t,t w ) which according to the above discusssion is a measure 
of inhomogeneities in the relaxation process under consideration. 



B. Trap models 

As already mentioned above, Bouchaud and Dean[2I] derived eq.tJH} for a trap model, 
although with £(t,t w ) = 0. In order to understand this, one has to consider the quantities 
of interest in the trap model in more detail. Instead of the correlation function given in 
eq.(jlj), Bouchaud and Dean consider a slightly different function: 

C(t,t w ) = q EA U(t,t w ) 

where U(t,t w ) denotes the propability that the system has not left the initial trap during 
t + t w and qEA is the Edwards-Anderson order parameter EBj. Note that U(t, t w ) can be 
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interpreted as a structure factor for large wave- vectors, because every escape out of a given 
trap gives rise to a decorrelation and thus to a decay of H(t, t w ). If one now identifies the 
unspecified states k of the last section with the traps and neglects any correlation between 
the values M k and the trap energies, only the G kk (t) are relevant G5j ■ Assuming that the 
M k are distributed among the traps according to some fixed distribution p(M), one finds 
for the correlation function 

c(t, t w ) = n(t, t w ) = (M 2 ) ]T G kk (t - t w ) Pk (t w ) 

k 

with qEA = (M 2 ) and accordingly 

i(t, t w ) = (M) 2 ]T Gkkit - t w ) [W klVl {t w ) - W lkPk (t w )} 

k,l 

From this last expression it is evident that for a symmetric distribution p(—M) = p(M) 
£(t,t w ) vanishes, £(t,t w ) =0, and the relation between the response and the correlation 
reads as 

R(t,t w ) = (3 [ad tw C{t,t w ) - (1 - a)d t C(t,t w ) 

For a = 0, i.e. a dependence only on the initial trap, this expression reduces to the one 
given by Bertin and Bouchaud|27j. Very recently, Sollich[23| has shown that the above 
expression holds for any trap model independent of the form of the transition rates. It 
will be shown later that £(t,t w ) = is also found for other models in which any transition 
yields a complete decorrelation of the stochastic variable of interest. 



C. Composite Markov processes 

In this section I will consider the following class of models. In order to relate the quantities 
M k to the values of k, I now treat the states as configurations of the system, characterized by 
their (free) energies e k and consider the composite Markov process (CMP) {e(t).m(t)}\17\. 
Typically, one is interested in the dynamics of the process m(t) only. For example, when 
considering the relaxation properties in a (free) energy landscape determined by some 
model, often a random assignment of the m k to the states k is used. The stochastic process 
m(t) is obtained from the CMP {e(t),m(t)} via the definition of marginal probabilities, 
i.e. by integrating out e(t)[T7]. Note, that usually the process m(t) defined this way is no 
longer a Markov process. Generally, m(t) can be related to the quantities M k in various 
ways. For instance, in the magnetic models to be treated in the following section one 
identifies the M k with 'magnetizations' m k . Another example is given by the case of the 
reorientation of a molecule in a liquid and the application of an electric field. In this case 
one chooses M k =pcos(Q k ) with fl k denoting the orientation of the molecule and p the 
value of the static dipole moment. Assuming that Q is a random variable, one considers 
the CMP {e(t), Q(t)}. Similarly, one can treat the translational motion of tagged particles 
by identifying M k with e* qrfe , where denotes the position and now the CMP {e(t),r(t)} 
is considered. 
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The CMP {e(t),m(t)} is completely specified by the transition rates W(e k ,m k \ei,mi). 
Of course, there are many ways to choose these rates and in the present paper I will make 
the specific assumption: 

W(e k ,m k \ei,mi) = W(m k \mi) +K W A( W )(m fc |mj) (12) 

Here, the W(m k \mi) denote transition rates for a mi—>m k transition and the K k i those for 
a q — > e k transition. Finally, the A( k i)( m k\ m i) determine what kind of change takes place 
with the m k in case of a e\ ^e k transition. The choice made in eq. ()12|) is general enough 
to cover most situations of interest. Furthermore, it will be assumed throughout that 
no correlation between the e k and the m k exists initially and therefore the corresponding 
probabilities factorize: 

p({e k ,m k };t = 0) = p(e k ; 0)p e k g (m k ) = p° k p e k q (m k ) (13) 

This means that additionally only the states k are affected by the initial non-equilibrium 
situation. This choice is sufficient for all models considered in the present paper. 

The two classes of models that will be considered in the present paper are completely 
determined by the choice of the W(m k \mi) and the A^(m k \mi) . Random magnetic models 
can be obtained from eq. ()12j) with the choice 

W{m k \mi)=ti and A( kl )(m k \mi) = a k (m k ) (14) 

where (j{m k ) denotes the a priori distribution of the m k in state k, usually assumed to 
be independent of k, a k {m k ) = u(m k ) Wk. Note that A^ k i){m k \mi) = a k {m k ) means that 
this quantity solely depends on the destination state. The implications of eq.(JT31) will be 
worked out in the next section. 

Another class of models concerns the models for dipole reorientation and translational 
noted already above. Here, one considers the CMP {e(t),a(t)}, with a(t) = Q(t) or a(t) = 
r(t). In this example, eq. (fT2|) reads as: 

W(e k , a\e h a) = W k (a\a')5 K i + K k iA(a\a!) (15) 

where I have already neglected a possible dependence of A on k, I. Two classes of models can 
be distinguished in the present context. One class is given by the so-called environmental 
fluctuation models and is defined by the choice of the on-site relaxation rates W k (a\a') 
and of A(a\a'). Usually, the latter quantity is either set to unity (no particle movement 
associated with an exchange) or one chooses A(a\a') =p eq (a) (complete randomization). In 
case of rotational motion the latter choice means that with any exchange process a random 
reorientation takes place [2.81 ESI ED]- 

The 'energy landscape model' is defined by assuming that there is no on-site relaxation 
at all, W k (a\a') = 0, and the geometry of the molecular motion is determined by A(a|a')[SI]. 
Note that in this model the process e(t) really defines the dynamics and a(t) can be viewed 
as some kind of slave process. In its applications to the dynamics of supercooled liquids 



7 



the assumption of angular jumps of a finite width Af2|32], i.e. A(f2|fi') = S(Q — [Q' + AO]), 
has proven successful [33 , but also random rotational jumps, A(Q\Q') =p eq (Q), have been 
used for general considerations |3*T] . In case of translational motion of tagged particles, an 
isotropic model with A(r|r') ocS(r — [r 1 + SR]) with a jump length 5R has been used to give 
a simple explanation of the apparent translational enhancement in supercooled liquids [3Uj. 
Before considering these kinds of models further, I will treat the random magnetic models 
in the following section. 

III. Random magnetic models 

In the usual treatment of magnetic models one simply assumes that a random magneti- 
zation m k is assigned to the states 'fc'|T[Jl El- Here, I will treat the CMP {e(t),m(t)} as 
introduced in the preceeding section with the transition rates determined by eq.lflljl. i.e. 

W(e k ,m k \ei,mi) = a k (m k )K M (16) 

The initial condition, eq. (fK?j) . for these models reads as p(e k ,m k ;0) = p k <Jk(rn k ). Start- 
ing from the ME, eq.(j2J), where now J2k nas to be read as J2kfdm k , one finds that the 
ansatz G({e k , m k }, t\{e h mi}, 0) = cr k (m k )G k i(t) with G M (t) = G(e k , t\e h 0) solves the ME. 
The same applies for the time dependent probabilities, P({e k ,m k },t) = a k {m k )P k {t). In 
calculating the correlation, the response and t w ) some care has to be taken in per- 
forming multiple summations. For instance, in an abbreviated form one has £(t,t w ) = 
J2k,i,n Jdm k fdmi Jdmiak(m k )ai(mi)a n (m n )m k miA k j^ n after all substitutions have been 
performed. In such an expression, one has to treat the terms k = I and k ^ I sepa- 
rately and then perform the integrations. This way all quantities of interest can be ob- 
tained from eqns. ()4p8p9pi0|) . In the most general case the results depend on the moments 
(m k ) — J dm k crk(fn k )m k . For example, the correlation is found to be given by: 

C(t,t w ) = ^2(Aml)G kk (t -t w )p k (t w ) +^2{m k )(mi)G k i(t - t w )pi(t w ) (17) 

k k,l 

with (Am|) = ((m k — (m k )) 2 ). If even distribution functions, a k {m k ) =cr k (—m k ), for which 
(m k ) = are chosen, the second term in the above expression for the correlation function 
vanishes. For £(t,t w ), one finds in this case: 

£(t, t w ) = J2(ml)G M (t - t w ) [ni k p k (t w ) - K k ipi(t w )) (18) 

k,l 

These expressions can be analyzed for arbitrary random magnetic models. It has, however, 
to be pointed out that the simple choice, eq. f)lfij1 is meaningful only, if the magnetizations 
are randomly distributed among the energies, because this cannot account for any corre- 
lations between the m k and the e k . In the following, I will illustrate the general formulae 
for the case of a simple model for which the master equation can be solved analytically. 
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Example: kinetic random energy model 

As an example, I will treat the special case of a kinetic random energy model (REM) 
with a special choice of the transition rates K k i- I will follow the treatment of Koper and 
Hilhorst [33] , who considered various forms of the transition rates. For the purpose of the 
present paper the simplest choice 

K kl = K exp {-(3e k ) (19) 

is sufficient, where k is a rate constant, to be set to unity in the following, and (3 = T~ 1 . 
This means that the k^i only depend on the destination state of the transition via the 
Boltzmann factors = exp (—/?€&). The corresponding ME can easily be solved giving 
G ik (t) = Z((3)- l Bi + \5 ik - Z((3)- l Bi]e^(-Z((3)t) with the partition function Z{(3) = 
J2k^k and 5ik denotes the Kronecker symbol. The equilibrium populations are given by 
p? = Z(p)~ 1 B i . 

In principle one can choose the random variable to depend on the e&, a possible choice 
being (ml) = B^ for arbitrary n, cf. the discussion given on this point by Fielding and 
Sollich 3 . Without showing explicit results here for the case n^O, it is appropriate to 
mention that also for the simple model considered in this section, the FDT-violations 
depend on the value of n, meaning that the m k chosen this way do not represent neutral 
variables PI In the following, I will choose n = 0, (m k ) = and (rn|) = 1. In addition, 
the initial populations p° will be chosen as p® = N~ 1 , appropriate for a quench from high 
temperatures in the beginning of the experimental protocol. For this case one finds, using 
the abbreviation Z{(3) = Z(2(3)/Z(f3) 2 : 

(t, t w ) = tfCf> (t-t w ) + AC< B > (t, t w ) 
Cg) (r) =Z{0) + [\- Z{0)\ exp (-Z09)t) (20) 
AC {B) (t,t w ) = Z((3) [exp (-Z(P)t) - exp (-Z(f3)t w )} 

and 

For the response one has 
R {B) (t,t w ) 

In these expressions, the explicit dependence on the Boltzmann factors Bi has been denoted 
by the superscript B. Interestingly, the expression for the response is independent of a. 
It should, however, be noted that in the more general case of a variable (m\) =B^, n^O, 
R( B )(t,t w ) is found to depend on a. 

Below the transition temperature T c = J / (2-y/log 2) the random energies are exponen- 
tially distributed [34 j, where J is related to the variance of the Gaussian distribution of the 
REM[35j. The corresponding distribution of Boltzmann factors B^ is given by 

p (B) = -xB- 1 - x for <B<oo (23) 



(21) 



dC%\t-U 
dt 



Z(P)Z(P)exp(-Z(f3)t) 



(22) 
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and p(B) =0 otherwise. Here, x = T/T c , v = l/(2-y/7r log 2) and N denotes the number of 
random energies. All quantities of interest are obtained as averages over the distribution 
P{B), i.e. 



F(t) = (F^(t)} 



B 



(*) 



(24) 



In order to perform the averages of the quantities given in eqns. (j2QI21l22J) . the only integrals 
required are: 



$(t) = ( e - z(/3) *> = exp 
V(t) = (Z(f3)e- zwt ) 



-vt x ) 



:i-x) 



$(t) - v 1/x T(l - l/x;vt x ) 



(25) 



where v = vT(l — x) with the Gamma function Y(a). Furthermore, T(a,b) denotes the 
incomplete Gamma function. The calculation of the averages follows the lines of ref.j 
and poses no problem. 
In particluar, one finds: 



C eq {r) 



Mr) 



1 - x + $(t) 
R(t,t w ) = -pd t [C eq (t-t w )+V(t)) 



(26) 



and 



at,t w ) = d tw y(t w )-d t y(t) (27) 

Koper and Hilhorst |H4| treated C eq (r) in their investigation. Using the asymptotic behavior 
of r(a, b) it can be shown that this function behaves as C eq {r) ~exp (— vt x ) at long times. 
For the integrated response xifiiw) — ft dt'R(t,t') one finds: 



X (t, t w ) = (3 [C eq (0) - C eq {t - t w ) - (t - t w )d t V{t)) 



(2f 



where the first two terms correspond to the equilibrium response, Xeq( T ) = ~ C eq (r)], 
because C eq (0) = 1 for n = 0. The limiting behavior of the integrated response can be 
summarized as follows: 



t 



oo 
-> 
-> 



X(t w + r, t w ) 

Xifw ^w) 
X\Pw tw) 



Xeq(r) 

(3x [1 - $(r)] Vr 



Vr 
Vt,„ 



(29) 







This behavior is illustrated in Fig.l, which shows a typical FDT-plot|3j, /3~ 1 x(t w + T,t w ) 
versus C(t w + r,t w ), for x = 0.3. From this plot and eq. (J2HJ it is evident that the limiting 
fluctuation-dissipation ratioP for this model is given by = T~ l . The dotted lines in 
Fig.l represent a slope of (—1) (ordinary FDT) and a slope of (—x), the limiting behavior 
for long times r. Therefore, in this simple model, the effective temperature coincides with 
the transition temperature T c . 
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The function [— £(t w +T, t w )} is plotted versus r in Fig. 2 for various values of the waiting 
time t w for x = 0.3. The behavior for other values of x is very similar. For very short r 
the function £ starts from zero and approaches a plateau value determined by [dt w ^{t w )} 
for long times r, cf. eq.(J27j). Only in the limit of long waiting one has £(t w + r,t w ) =0. 
The lower panel of Fig. 2 shows a logarithmic plot of (— £), which demonstrates the linear 
behavior of £ and also that the cross-over to the plateau value happens around the waiting 
time, i.e. for r ~ t w . For not too long waiting times, the approximate expression 



holds with high accuracy. This expression is shown in the lower panel of Fig. 2 as the thin 
dotted lines. For t w = 10~ 5 , it cannot be distinguished from the exact expression, eq. (|2*7jl . 

From these considerations it becomes clear that even for this simple example the as- 
sumption that the distribution function determining the relaxation is equilibrated with 
respect to the e*, is meaningless. On the contrary, the system gets stuck in a distribu- 
tion characterized solely by the transition temperature T c . In order to further discuss the 
properties of the function £(£, t w ), in the next section I consider a different class of models. 

IV. Models for tagged particle motion 

Models exhibiting truely non-equilbrium dynamics are often considered as models for glassy 
systems. However, when one is concerned with supercooled liquids, the usual situation 
is one in which the system can reach metastable equilibrium for long times. The glass 
transition temperature is merely a convenient way of classifying the time scale of the pri- 
mary relaxation |36j. Nevertheless, investigating the out-of-equilibrium dynamics of simple 
models for molecular motion in supercooled liquids is capable of providing additional in- 
formation about the mechanisms of the dynamics. 

In this section, I consider models for tagged particle motion that have been used suc- 
cessfully to describe the translational and reorientational dynamics in supercooled liquids 
on a phenomeno logical level [2ffi ETT] 153]. 

In case of the reorientational motion, the correlation function C(t,t w ), eq.(jl|), is ob- 
tained by identifying M(t) with orientation dependent static dipole moments fJ.(t). The 
time dependence has its origin in the temporal fluctuations of the direction of fi relative 
to a space fixed axis. It should be pointed out that C(t,t w ) is directly related to the so 
called rotational correlation function gi(t,t w ), where the general definition is 



Here, the angle 6 is the polar angle of the axially symmetric principal axes system in the 
laboratory fixed coordinate system and Pl{x) denotes the Legendre polynomial of rank L. 
For instance, in dielectric spectroscopy 6 is the angle between the applied field and // and 
gi(t,t w ) is measured. In NMR or light scattering g2(t,t w ) typically is observed. 




(30) 




(31) 
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If one is interested in translational motions of tagged particles, one choice is to identify 
M(t) with e iqr (*), with q denoting the wave- vector and r(i) is the position of the particle. 
This means, one has a close relation between C(t,t w ) and the incoherent intermediate 
scattering function |37]: 

C^{t,t w ) = ( e *v(*) e -*irf*»)) (32) 

where I have already assumed that the system is isotropic, q = |q|. 

In the following I consider the CMP {e(t),a(t)}, where the a(t) either denotes the 
orientation Q(t) or the position r(t). 

Energy landscape models 

I start the discussion with the CMP {e(t),a(t)} in the context of the energy landscape 
model introduced in ref.|23- ^ n this model, one associates a (free) energy with each 
basin or metabasin of the potential energy landscape of the system^ I38j. In the model 
of a random first order transition put forward by Wolynes et al., see e.g. ref.jHH], the 
would be interpreted as the free energies of the metastable glassy states, characterized by 
an analogue of the Edwards- Anderson order parameter. The slow primary relaxation in 
such scenarios is determined by the transition among the various states. In the energy 
landscape model a tagged particle is allowed to reorient about an average jump angle 
G and also to hop a mean distance SR whenever a k — > I transition takes place pH] EJ3]. 
As already mentioned in Sect.IIC, the simplest choice for the transition rates given in 
eq.(|T5|l is Wk(a\a') = and A(a\a') oc 6 (a — [a' + Aa\) with Aa denoting 6 or 5R. In the 
context of reorientational motions, sometimes it is sufficient to consider rotational random 
jumps, which just means that any k — > I transition randomizes the molecular orientations 
completely. In this case one chooses A(Q\Q') = l/(8ir 2 ). 

In Appendix A the solution of the ME for this model is explained and the expressions for 
the quantities relevant in the present context are derived. The correlation function C(t, t w ) 
is given in eq. (jA.6|) and £(t,t w ) in eq. (jA.7|) . From the latter expression it is obvious that £ 
is closely related to the non-equilibrium part of the correlation function, AC. The response 
then follows from the general expression, eq.pjl. For £(t,t w ) it is found that it is given by 

£ {1) (t,t w ) oc cos(0) 

in case of reorientations with an jump angle 0. If random reorientations are assumed 
instead, one has t w ) = 0. For translational motion in an isotropic model one finds 

£ (Q) (Mjocj (Q) 

where jo(x) is a Bessel function and Q = (q ■ 5R) with the wave- vector q. The limiting 
behavior is found to be given by {p = Q or L): 



t n 



oo : 




+ T, t w ) - 


^0 


Vr 


oo : 




+ r, t w ) - 


^0 


vt w 


> : 


£ {p) (t w 


+ t, t w ) - 


-£ (P) M) 


Vr 


> : 




+ r, t w ) - 


^0 





(33) 
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It is evident that in contrast to the simple kinetic REM considered in the previous section, 
£(t w + T,t w ) vanishes in the limit of long r. The reason for this difference is that in the 
present model the system reaches equilibrium in this limit. 

Although the general expressions for £ are not very instructive, two limits can be an- 
alyzed in simple terms and the relevant formulae are given in Appendix A. One limit is 
given by either random rotations or large wave-vectors. In both cases, one finds that 
£ vanishes. This means that whenever a transition gives rise to a randomization of the 
relevant variable, £ is expected to vanish. This is very similar to the case ot the trap 
models discussed in Sect.IIB. Therefore, in this situation, the correlation function C(t,t w ) 
is essentially the same as the function U(t,t w ), giving the probability for staying in the 
same orientation/position (as opposed to the same trap). This can also be seen from the 
direct comparison of the expression for C^(t,t w ) and C^\t,t w ) given in eq. (jA.9|) with 
the one given for U(t, t w ) in Sect.IIB, because only G/^t) =e~ Kfct enters these expressions, 
cf. eq. (jA.8|) . The correlation functions in this limit are independent of L or Q because 
every transition gives rise to a complete decorrelation. Thus, there is a close link between 
the simple model considered here and trap models concerning the vanishing of the function 

t w ). It should, however, be noticed that this correspondence relies only on the 'geomet- 
ric' aspects and has nothing to do with the aging dynamics in the models considered. The 
important point regarding t w ) is the fact that due to the randomization of the relevant 
stochastic variable in each transition the distribution of orientations/positions is homoge- 
neous. A very important difference between the model for translational model and the trap 
model, however, is that the wave- vector is a quantity that can be adjusted experimentally. 
Therefore, the limit of vanishing £(®> can be reached continuously in principle. 

The other linit that can be investigated in detail is that of very small rotational jump 
angles (rotational diffusion) or small wave-vectors q. In this case, one finds that the 
equilibrium correlation functions decay approximately exponentially in time with diffusion 
coefficients related to the averaged transition rate (k) =J2k,i K kiPk , cf. eq. (jA.llJ) 3.3*]. In 
this limit, the prefactor of £, cos (0) or jo(Q), are approximately equal to unity. 

The general expressions for the correlation function and the function £ are illustrated 
in Appendix A by a simple two-state model. As is to be expected, AC and also £ vanish 
for equal transition rates, K\2 = K21, implying pl q = p e 2 q ' ■ Also the transient nature of AC 
and £ is obvious from eq. (|A.13|) and eq. (|A.14|) . 

Environmental fluctuation models for molecular reorientations 

In order to show that the results obtained in the framework of the energy landscape model 
qualitatively are the same also for other models, I will now discuss the exchange models 
mentioned earlier. For simplicity it is asumed that there are two 'states', one in which 
molecular reorientation is fast, denoted by '/' i n the following, and one in which molecules 
reorient slowly, denoted by Vj2H|- The generalization to an arbitrary number of states 
poses no problem and has been used to describe several features of the reorient at ional 
relaxation in supercooled liquids, including dielectric relaxation (HO) and four dimensional 
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NMR experiments |29| I4*U] . If the exhange rates, K s f and Kf s , vanish, one has a static 
heterogeneous dynamics in the sense that slow molecules remain slow forever. The solution 
of the corresponding ME and the relevant results are given in Appendix B[2E], where it is 
shown that there are only slight differences to the case of the energy landscape model. 

For the case of no changes in molecular orientation due to k — ► I transitions, orien- 
tational relaxation takes place solely due to on-site reorientations. From the expression 
for AC(t,t w ), eq. ()B.2j) . it becomes evident that this function vanishes, if the reorientation 
rates in the two environments, Tf and T s , are equal, i.e. AC(t, t w ) = t w ) = for r s = T/. 
The same holds if one has K s f = Kf s . Furthermore, it is evident that in this model the only 
relaxation during the waiting time is given by exchange. For vanishing exchange rates £ 
vanishes, but AC remains at its initial value for all waiting times. The latter fact is easily 
understood because in this case there is no way for equilibration of the V° s /f- As in the 
energy landscape model, £(t w + r, t w ) vanishes in the limits of small and large r and has a 
maximum at intermediate r determined by the effective decay rates. 

Sofar, I have considered the two-state model in which an exchange process has no effect 
on the molecular orientation. However, the other extreme scenario, in which the molecular 
orientation randomizes completely in case of a k — > I transition, can easily be treated in the 
same way. It turns out that the only relevant difference to the model considered above is 
that now £(t, t w ) = 0, cf. eq. (jB.4|) . In this case the orientational distribution is homogeneous 
at each instant of time and the argument for a vanishing £ given in Sect. IIA applies. 

V. Conclusions 

I have considered the fluctuation-dissipation relation for a general class of stochastic models 
obeying a master equation. I have chosen transition rates which in the presence of a field 
are disturbed in a multiplicative way. Only cases for which detailed balance is obeyed with 
and without field were considered. 

The main result of the present paper is given in eq.tjHJ) and shows that in the general 
case the response cannot be related to time-derivatives of the correlation function alone. 
Instead, a function £(t,t w ) occurs additionally in the expression for the response, which 
cannot be related to any physical quantity in a simple manner. An exception to this finding 
is provided by trap models, where it has been shown the t w ) vanishes for even distribu- 
tions. Hoffmann and Sibani[T8t IT9] were the first who considered the function £(t,t w ) for 
the special case of a = 1 in the general expression for the fluctuation dissipation relation, 
eq.fJHJ). However, they did not further discuss the meaning of this function apart from the 
fact that it vanishes if the relaxation is determined by probability distributions that are 
homogeneous with respect to the states of the stochastic process under consideration. 

In the present paper I have shown that this assumption is likely to be wrong even for 
very simple minded models. As one example I considered a kinetic random energy model 
with an extremely simple choice for the transition rates in Sect. III. The limiting effective 
temperature was found to coincide with the transition temperature of the REM, T eff = T c . 
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It was found that even for this model £(t w + r,t w ) ^ for t w < oo. Due to the fact that 
the system never reaches equilibrium in this model, one also has a non-vanishing long time 
limit £(oo, t w ). 

Furthermore, I have shown that for some models in which any transition leads to a 
complete loss of correlation of the stochastic variable under study, one finds that the 
function £(t,t w ) vanishes. This holds for trap models and also for models of molecular 
reorientations, provided these reorientations proceed via random rotational jumps. If the 
the energy landscape model is applied to translational motions of tagged particles, one 
has (,( Q \t,t w ) oc jo(Q), Q = q-SR, if the intermediate scattering function C^(t,t w ) is 
measured. If Q ^> 1, one therefore also finds that vanishes. In this example, the 
fact that £ vanishes thus relies on an experimentally adjustable parameter and does not 
reflect a property of the transition rates. It is well known, that a stimulated NMR echo 
experiment is the 'rotational analogue' of the intermediate scattering function|41j. If this 
experiment is performed in a way that is analogous to C^^'it, t w ), every single rotational 
jump leads to a decorrelation and the random rotational jump model applies. Also in case 
of the trap model any transition out of a trap leads to a loss of correlation. Thus, also 
n(t, t w ) corresponds to C^ Q2>1 ^(t, t w ). Therefore, it is tempting to speculate that whenever 
an analogue of C^'(t,t w ) in the large Q-limit is considered for C(t,t w ) in the sense that 
every transition gives rise to a complete decorrelation, the function £^(t,t w ) is expected 
to vanish. This also means that whenever every transition among the states of a system 
obeying a ME randomizes the stochastic variable under consideration, I expect £(£, t w ) =0. 

For all models considered in the present paper, it turns out that £(t,t w ) is related to 
the non-equilibrium part of the correlation function, AC(t,t w ). Furthermore, the main 
difference between the kinetic REM and the models considered for tagged particle motion 
consists in the fact that in the latter models the system reaches equilibrium for long times. 
It is found that £(t w + r,t w ) (if it is non-vanishing) for these models is of a transient 
nature and vanishes for short and long r, in contrast to the case of the kinetic REM, where 
lim^oo + r, t w )^0. 

At the moment it is not clear, how the function £(t,t w ) behaves for other types of 
models. A study of kinetic Ising chains is under way. Also, it would be interesting to 
investigate the behavior of £(t,t w ) for varying initial conditions. In the present paper, 
eq.([T3]) was always assumed to hold. 

To conclude, I have shown that for some general class of models the fluctuation dissipa- 
tion relation is determined by time-derivatives of the correlation function and an additional 
function £(t,t w ). As noted earlier [TH]. £(t,t w ) is a measure for deviations from homogene- 
ity in the distribution functions determining the relaxation. This function is expected to 
vanish in situations where transitions among the states considered lead to a randomization 
of the relevant stochastic variable, meaning that in such a case the relevant distributions 
are homogeneous. 
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Appendix A: Energy landscape model 

In this Appendix I will give the relevant formulae that are needed for the solution of the 
ME for the CMP {e(£), a(t)}, where a(t) has to be identified with the orientation Q or the 
position r, depending on the situation. Because the solution of the ME proceeds in the 
same way for both cases, the treatment can be formulated quite generally. The transition 
rates are given by 

W(e k ,a\ei,a!) = K k iA(a\a') (A.l) 

which follows from eq. lfTK]) for Wk(a\a') = 0. Again, I neglected possible dependencies of 
A(a\a') oc 5 (a — (a' + Aa)) on the initial and final states of the transition. To solve the 
corresponding ME, eq.(J2J), for this case, one defines the matrix of the eigenvectors of A(a\a'), 
U(a,p) corresponding to the eigenvalue A(p), i.e. A(p) = Jda Jda'U~ 1 (a,p)A(a\a')U(a,p). 
Next, the conditional probability is expanded in terms of these eigenvectors, 

G({e k ,a},t\{e h a },0) = j "dpU(a,p)G% \t)U-\a Q ,p) (A.2) 

The Greens functions (t) are then found from the solution of 

G&\t) = -K„G%\t) + A(p) £ KknG<j$(!t) (A.3) 

n 

where the sum rule K k :=J2tek K ik was used for the diagonal element. The initial populations 
are chosen according to eq. (|13j) . which in the present context reads as: 

p({e k ,a};0)=P e9 (*)Ph (AA) 

If one considers isotropic rotational motions with a fixed mean jump angle G, one has 
to identify a with Q. In this case one has p eq (a) = ^y, U(a,p) = J^^D<£l(Q), A(p) = 
Pl(cos (0)) and the integration over p is now a sum over L. 

In case of translational jumps onto all postions of a sphere with a radius 5R (jump 
length), the corresponding substitutions in the general formulae are p eq (a) = V^ 1 with V 
denoting the volume. Furthermore, one has U(a,p) = ^e* qr and A(p) = jo(q5R), where 
jo(x) denotes the Bessel function of zeroth order. 

In eq. (jA.3|) one has to identify p = L for rotations and p = Q with Q = (q-SR) for 
translations and thus eq. (|A.3|) reads as: 

Rotation : G^(t) = -KyG$\t) + P L (cos (6)) ]T K kn G ( ^(t) 

n 

Translation : G^(t) = -« k G^ Q) (t) + j (Q) £ n kn G^\t) (A.5) 
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Note that due to Po(x) = jo(0) = 1, there is one eigenvalue, p = 0, A(0) = 1. Therefore, for 
p = 0, from eqns. (jA.31 IA.5|) the original ME for the transitions in the energy landscape is 
obtained, i.e. (t) = Gki(t). This is in full accord with the idea underlying the model 
that the variable a is allowed to change only in case of a k — > I transition among the states 
of the system p3: . 

Once the G$(t) are obtained from the (numerical) solution of the equations (IA.5J1 . all 
quantities of interest can be calculated and one finds, cf. eq.©: 

^M = EGlf(r) P r ; AC^(t w + r,t w ) = j:G^(r)[pi(t w )-pt q ] (A.6) 

k » L & jit 

and 

Z {p) (t w + r,t w ) = A(p)^- (AC^(t w + r,t w )) (A.7) 

In these expressions the Pk(t w ) = Y^iGki{t w )p < j l . This is easy to understand from the fact 
that according to the chosen initial conditions, eq. (|A.4J) . the variable a was in equilibrium 
in the beginning, t w = 0. Furthermore, in case of reorient at ional motions I have omitted 
a trivial factor |/i 2 . Note that for dipole reorientations, one is concerned with C^(t,t w ) 
and ^'(tjtw), i.e. one always has to use L = l in the present context. 

If in case of a = Q it is asssumed that the reorientations proceed via random jumps, one 
has to replace the factor Pi, (cos (0)) in eq. (|A.5J) by 5l,o, meaning that for L > only the 
diagonal element is retained. 

In general, the behavior of the function ^\t,t w ) and all other quantities is obtained 
from a numerical solution of the equations for the Greens functions, eq. (|A.5)) [ 3*T) I33j. For 
this purpose, the transition rates km have to be chosen in a prescribed way. It is however, 
possible to give some general results valid in specific situations. Consider the case of 
rotational random jumps, for which one has A(p > 0) = 0. Similarly, if the incoherent 
scattering function is observed for large wave- vectors, Q > 1 (or g > (5R) _1 ), one has 
jo(Q) — 0. In both situations it is clear that £^(t, t w ) vanishes. In the other extreme, 
namely very small wave- vectors or rotational diffusion, one has A(p) = l—f]p with t]q = Q 2 /6 
and tjl = L(L + l)(0/2) 2 . In both of these limiting situations, one can treat the problem 
in perturbation theory, as outlined below. 

A(p)^0: 

In this limit, the general expression for the Greens function, eq. (jA.3|) . approximately reads 
as Gffity) = —KkG^i(t) and therefore one has {nk'-=Hi^k K ik) 

G%\t)~6 u e-"* (A.8) 

Using Apk(t w ) =pk(t w ) — p e k q , one easily finds from eq. (|A.6|) : 

C%\t) -J2 e ~ KkT Pk and AC^(t w + r,t w ) ~ 5>~^A Pk (t w ) (A.9) 

k k 

From ea. flATJ) it is obvious that ^ p \t,t w )^0 for A(p)->0. 
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A(p) = l-77 P , 77 P <1: 

In this case one is in the diffusive limit and one has t)q = Q 2 /6 and t]l = L{L + l)(@/2) 2 , as 
already noted above. Here, one can use the ME for p = as the unperturbed problem and 
consider perturbation theory with a perturbation given by matrix elements Vm = —r) p Kki- 
This can be seen directly from the ME, eq. (IA.3|) and using A(p) = 1 — 1] P therein, cf. ref.[33 . 
The unperturbed Greens functions thus coincide with the Guit) and up to linear order in 
rjp one has: 

G$\t) ~ G kl (t) - Vp J2 fdsG km (t - s)K mn G nl (s) (A.10) 

m,n 

A straightforward calculation yields: 

C ip) (t w + r,t w ) ~l-rj p ^(k)t + jf dsJ2 K k[Pk(s + t w ) -p*j?]J 
where the mean relaxation rate is defined by 



( K ) = J2 KkP e k 
k 

Therefore, up to linear order in rj p one finds 

C<p\t) ~ e ~^ T (A.ll) 

and 

AC®{t w + r, t w ) ~ - Vp [ T ds J2 «* [Pk(s + t w ) - p e k q ] e-^ T (A.12) 

Jo k 

Therefore, in the limit of small wave-vectors the model predicts an exponentially de- 
caying intermediate incoherent scattering function, C^(r) = S(q, r) ~exp (—q 2 D T r), with 
an apparent diffusion coefficient = (8R) 2 (k) /6|33|. In case of the rotational diffusion 
of molecules one accordingly has C^(r) ~ exp (— L(L + \)D r t) with D R = (6/2) 2 (ft). 
Remember, that in the present discussion only L = l is relevant. Again, £(t,t w ) is given by 



Two-state model 

In order to illustrate the results obtained above for the energy landscape model, I now 
consider only two states. This is sufficient in order to discuss the qualitative features of 
the general model. The elements of the Greens function matrix are explicitly given by 
G n (t) =pt 9 +p?e-« G 22 (t) =pf +pt q e-*\ G l2 (t)=p eq {l - e -«) and G 21 (t) =p\\\ - 
e~ Kt ). Here, I defined R = k± 2 + k 2 i- Furthermore, one has Api(t w ) = [p\p e 2 — p 2 pl 9 ] e~ Ktw 
and Ap 2 (t w ) = — Api(t w ). With these quantities one finds for the two limiting situations 
discussed above: 

A(p) -> : C$ (r) ~ e~ K21T pt q + e~ K12T p 2 q 

AC {p \t w + r, t w ) ~ \p° lP e 2 q - p° 2 p eq ] \e- K2lT - e~ Kl2T ] e~ Rtw (A. 13) 
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One therefore finds that vanishes in the limit of small and large measuring times r as 
well as for long waiting times t w . In addition has a maximum as a function of r for 

T max = \ll («12/«2l)/(«12 - «2l). 



A(p) 



*7p 



»7 p (k}t 

K 12 — K 2l 

R 



(A.14) 



eg e« 

P1P2 - P2P1 



l-e 



Here, the situation regarding £W is analogous, the difference being that the maximum of 
^ p) (t w + r, t w ) now occurs at T max = R' 1 In (l + -^x) . 



Appendix B: Environmental fluctuation model 

Such models are defined by the transition rates given in eq.(|15|). where now the on- 
site reorientations are modelled via finite Wfc(f2|f2'). For simplicity, it is assumed in 
the following, that the reorientations proceed via random rotational jumps, W^f^Q') = 
— r^8(p, — Q') + rfc/(87r 2 ). In addition, it has to be quantified what happens to the molec- 
ular orientation in case of k — > I transition. As in the original model of Beckert and 
Pfeifer j2Sj, it will be assumed that either no change at all, A(f2|f2') = 1, or a random rota- 
tion, A(fl|fi') = l/(87r 2 ), takes place. Proceeding exactly in the same way as in Appendix 
A, the Greens functions are obtained from: 



(B.1) 



Here, I defined c L = l, if A(fi|fi') = l and c L = 5 Lfl if A(fi|fi') = 1/(8tt 2 ) j22j. As in the case 
of the energy landscape model, the relevant quantities are given by eqns. (jA.6|) and (jA.7|) 
if again it is assumed that there is no correlation initially, p({tk, Q}; 0) = -^Vk ■ 

In the following, the explicit expressions will be given for the simple two-state model 
discussed in Sect. IV. Here, I consider the case cl = 1. As in the text, the states are 
denoted by / and s. The equilibrium populations of the states are given by p e s q and p eq , 
the reorientation rates are T s , Ff and the exchange rates are denoted by K s f and Kf s . The 
solution of eq. (jB.l|) poses no problem and one finds, again omitting a factor /i 2 /3: 



c<£{r) 

AC^(t w + r,t w ) 



2w 

1 
w 



z f e~ A f T + z s e 



-A,r 



(B.2) 



„0 eq 
PfPs 



eq 
PsPf 



-A 



/ _ 



2 t\>~t<u 



Using the abbreviations K = (l/2)(n s f+Kf s ), T = (l/2)(Tf+T s ), 7 = (l/2)(Tf—T s ) and A K - 



(l/2)(Kf s — K s f) one has for the effective decay rates Aj iS = f + R ± y (A K + 7) + Kf s K s f. 
Furthermore, I used w = y/'j 2 + 27 A K + R 2 and z// s = w =F R ± 7(j»j 9 — p e s q )- Finally, for the 
function £ one finds: 

& 1 \t,t w ) = -(2R)AC<- 1 \t,t w ) (B.3) 
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The most important thing that changes, if one assumes that a random reorientation is 
associated with every k— ±1 transition, cl — Sl^, in eq. (jB.l|) is that one now finds: 



£ (1) (MJ = for c L = 5 L>0 (B.4) 
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Figure captions 



Fig.l : f3~ 1 x(tw + T , t w ) versus C(t w + r, t w ), for x = 0.3 and for various waiting times t w . 
The dotted lines represent the slopes expected for ordinary FDT (slope: —1) and an 
effective temperature of T eS = T c (slope: —x). 

Fig. 2 : [— £(t w + r,t w )} versus r, for x = 0.3 for various waiting times t w . Upper panel: 
linear scale, lower panel: logarithmic scale. In the lower panel, the approximate 
expression, eq. pU|) . is shown as thin dotted lines. 
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G. Diezemann; Fig. 1 
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G. Diezemann; Fig 2 
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